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We study numerically, using a one-dimensional Heisenberg model, the spin-Peierls 
transition in the linear Cu^~^ spin- 1/2 chains in the inorganic compound CuGeOs 
which has been recently observed experimentally. We suggest that the magnetic 
susceptibility, the temperature dependence of the spin gap and the spin-Peierls tran- 
sition temperature of this material can be reasonably described by including nearest 
and next nearest neighbor antiferromagnetic interactions along the chain. We esti- 
mate that the nearest neighbor exchange parameter J is approximately 160 K, and 
that the next nearest neighbor exchange parameter is approximately 0.36 J. 
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The purpose of the present study is to describe the spin-Peierls transition in the hnear 
Cu^"^ spin-1/2 chains in the inorganic compound CuGeOs which has been recently observed 
|1|J^J^,^]. The transition temperature Tc ~ 14 K has been inferred from the rapid drop 
of the magnetic susceptibihty towards zero, indicating the opening of an energy gap for 
singlet-triplet spin excitations. ||^,^ The existence of this transition was also confirmed by 
measurements of the heat capacity [0 which exhibits a sharp anomaly at Tc corresponding 
to a second order phase transition. We adopt a simple model Hamiltonian consistent of 
antiferromagnetic Heisenberg interactions along a chain. Since the one-dimensional (ID) 
spin system has no phase transition at finite temperatures, it is necessary to take into 
account the spin-lattice coupling in order to describe the spin-Peierls transition. We consider 
the spin-lattice coupling in the adiabatic approximation. The study is performed by exact 
diagonalization on finite chains. We suggest that in order to obtain a reasonable fit of the 
magnetic susceptibility a next nearest neighbor antiferromagnetic interaction along a chain 
should be included in the model. If the predictions resulting from our model are confirmed 
by additional experimental work, the Cu02 chains in this compound would be one of the 
few experimental realizations of a ID spin-1/2 Heisenberg antiferromagnet with competing 
interactions. The present study should also be considered as part of a current theoretical 
effort to understand magnetic properties, in particular the singlet-triplet spin gap, in low- 
dimensional spin systems such as Sr2Cu406 and (VO)2P207. 

The thermal properties of spin-Peierls transitions in spin-1/2 antiferromagnetic chains 
have not been studied, to our knowledge, from a microscopic and numerical point of view. 
An analytical study based on the transformation of Pauli spin operators to spinless fermion 
operators has shown that the homogeneous magnetic chain is unstable with respect to dimer- 
ization as the temperature decreases and a spin gap appears as a result of such dimerization. 
p!0| , p!T|] The description of experimental data was usually performed using the mean-field ap- 



proximation on the spinless fermion Hamiltonian. |jT2|] The effect of a magnetic field on the 
spin-Peierls transition temperature has also been analyzed by Cross and Fisher |]T3[ using 
a Luther-Peschel-type treatment of the spin correlation functions. The predictions of this 



theory were compared with experimental results for CuGeOs, 0] confirming the spin-Peierls 
nature of the observed transition. 

The ID microscopic Hamiltonian for the spin degrees of freedom that we consider is: 

Hs = JiE 

§21-1 ■ S2i + J2 / , S2i • S2i+1 

(1) 

i i 

where the index i runs over the lattice cells (i=l,...,N/2, N: number of sites) with periodic 
boundary conditions. We assume linear dependence of the exchange integrals on the atomic 
displacements u, so that: 

Ji = J (1 + 7u) (2) 
J2 = J(l-7u) 

where 7 is a constant. It is convenient to introduce the dimensionless quantity 6 = 7U. 

The spin- lattice interaction enters in this model only through Ji and J2. The underlying 
physical picture is the following. The spin chains in CuGeOs are oriented along the c 
direction (see for example Fig. 1 in Ref. According to experimental results the 

relevant lattice distortions are observed along the b axis, perpendicular to the chain direction. 
In this situation, the variables u in Eq. () correspond to oxygen displacements below the 
spin-Peierls transition temperature. @] It has been recently reported |^ that there is a 
comparable shift of the Cu ions along the c direction. The exchange constants Ji and J2, 
calculated following the path Cu-O-Cu, have in principle a complicated dependence on the 
displacement u. Then, the expressions for Ji and J2 given above should be considered as a 
first order approximation in u (or 6) of this complicated function. Besides, the application 
of Hamiltonian ([l|) to this material implies that we are neglecting the two-dimensional (2D) 
interchain exchange interactions which have a magnitude of approximately 10 % of the 
intrachain coupling, [^j 

The main difference between the Hamiltonian given by Eqs. (^ and (), and the one 
corresponding to the dimerized or alternating bond chain, ||I3|,|T^ is that in the former both 
Ji and J2 are temperature dependent because 6 = S{T). We calculate the temperature 
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dependence of these couplings by minimizing at each temperature the free energy T of 
the total Hamiltonian H = Hg + Hph with respect to 5 (adiabatic approximation). Then 
^mm(T) = ^(5oq(T)). The elastic term of the Hamiltonian is: 

Hph = -NKu^ = 1^52 (3) 

where K is the elastic constant. It is also customary to introduce the dimensionless spin- 
lattice coupling constant A = J7^K~^. This method of calculation was devised by Beni and 
Pincus fl^ for their study of the spin-Peierls transition in a spin chain with XY interactions. 
For a finite chain, and for a given set of parameters J and A, once we have determined the 
equilibrium displacement b^q at each temperature, any thermodynamical quantity can be 
computed in the dimerized region. 

The first stage of our study consisted in estimating the parameters J and g by fitting 
the experimental data for the susceptibility [|l| with the theoretical curve in the uniform 
or non dimerized (Ji = J2 = J), region. We computed the susceptibility by generating all 
energy levels Ej and their multiplicities dj in each sector of fixed total S^, using a Householder 
algorithm. The susceptibility was then obtained through its relation to the expected squared 



magnetization, summed over energy levels and total sectors; 

X(T) = g /Xb/? . „_flF,. • (4) 



This approach has the advantage that explicit eigenvectors are not required. 

The observed average g-factor is approximately 2.14 with a slight anisotropy along the a, 
b and c-axis. The experimental data has a broad maximum near 56 K (see Fig. [l|). A fitting 
of this data using the uniform Heisenberg model reproduces the position of this maximum 
of the susceptibility if the exchange constant J is chosen to be 88 K. [|T| However, as it can 
be seen in this reference the overall fitting is quite poor. A somewhat better fitting can be 



achieved by choosing J = 170 K. However, for T < 150 K the fitting is still quite poor. 
We have seen numerically that there are no satisfactory fitting of the experimental data in 
the region T > T^ax = 56 K using a nearest neighbor Heisenberg model. One of the simplest 



extensions of this model is to include in the spin part of the Hamiltonian a next nearest 
neighbor interactions term: 

H2. = J'ESrSH2 (5) 
j 

where the index j runs over the lattice sites (j=l,...,N). The possibility of an antiferro- 
magnetic (J' > 0) second-neighbor coupling through the Cu-O-O-Cu exchange path was 
suggested in Ref. 0. Another possibility is to consider the interchain coupling which would 
lead to 2D model. Taking into account the underlying physical picture discussed above, it is 
reasonable to assume that J' is independent of u at least in first order approximation. This 



Ji — J2 — J' model was studied at zero temperature by Shastry and Sutherland |]T9[ for a 
particular relation of the parameters where the ground state can be exactly calculated. 

We determined J by imposing that the maximum in the susceptibility is at 56 K as 
indicated in Ref. fl]] . Then, we determined the ratio = J'/ J in order to fit the maximum 
of the susceptibility Xmax = x(Tmax)- For this fitting we chose the susceptibility measured 
in a polycrystal sample |18| shown in Fig. |l]. In this figure we also reproduce the magnetic 
susceptibility measured on a single crystal along the a, b and c axis. Taking into account 
this dispersion of the data, we conclude that a reasonable fit is obtained with the following 
set of parameters: 

J = 160K 
a2 = 0.36 

To simplify the calculations, we have taken g equal to 2.00. This value of g, which is some- 
what smaller than the experimental value, should be considered as an effective value since 
we are neglecting the interchain couplings. Notice that the temperature region where we are 
fitting the available experimental data is still far from the asymptotic regime described by 
the Curie law. The effect of the neglected interchain coupling on the fitting parameters has 
been discussed in the literature (See e.g. Ref. [0). The results for the calculated suscepti- 
bility and the experimental data are shown in Fig. [l|. In both theoretical and experimental 



data a contribution from the orbital part of the susceptibihty, x"^'' = 1.5 10~^ emu/mole 
has been added. In this temperature region, T > T^ax, the finite size effects are neghgible 
and aheady for N = 12 the resuhs do not vary by taking larger clusters. For temperatures 
smaller than T^ax but in the uniform region, i.e. above the spin-Peierls transition temper- 
ature, there are strong finite size effects. The magnitude of these finite size effects can be 
seen in the inset of Fig. |l| where we show the susceptibility below 30 K for N = 8, 10, 12, 14 
and 16, and N = 9, 11, 13 and 15. 

A possible unwanted feature of the Heisenberg model with nearest and next nearest 
neighbor interactions is the presence of a temperature independent spin gap, i.e., in the 
absence of dimerization. This spin gap can be inferred from an exact solution at zero 



temperature found by Majumdar and Ghosh [|20| and was confirmed by subsequent numerical 



work. Recent studies on this problem indicate that for a2 > 0.25 there is a finite singlet- 



triplet gap. 1^ For a2 = 0.36, we have calculated the spin gap at zero temperature on finite 
lattices with N < 24 spins. To extrapolate to the bulk limit we adopted the form predicted 
by spin-wave theory, or alternatively, the essentially equivalent law. 



A(N) = + ^ (6) 

The resulting extrapolated value of the spin gap is approximately equal to 0.015 ± 0.005 in 
units of J, or 2.4 ± 0.8 K, much smaller than the smallest measured value for the CuGeOs. 
i 

The second stage in the calculation was the estimation of the coupling constant A. In 
order to do this estimation, we chose another piece of experimental data, the singlet-triplet 
spin gap at zero temperature. The fitting of this data is very convenient from the numerical 
point of view since at T = we can diagonalize larger lattices using the Lanczos algorithm. 
Using the spin part of the Hamiltonian given by Eqs. (|1|) and (^), we first determined the 
value of 6 that reproduces the experimental singlet-triplet spin gap which is ~ 2.15 meV 
(from Ref. [Q) or 0.153 in units of J = 160 K. Results for the spin gap for several values 
of 6 and for several sizes are shown in Fig. |^a. The extrapolation of the spin gap to the 



bulk limit for each value of S was also done using Eq. (|^). The extrapolated spin gaps 
as a function of 6 are shown in Fig. A quadratic interpolation gives the final result: 
4g(T = 0) = 0.014 ±0.001. 

Then, for each lattice size N = 8, 10, ...,22 and for 5 = 0.014, we computed the inverse 
of the coupling constant A(N)~^ that minimizes the ground state energy of the total Hamil- 
tonian, Hg + + Hph. The results are shown in Fig. |]c. Finally, by extrapolating A(N)~^ 
to the bulk limit we obtain = 21.3 ± 0.5 or A ~ 0.05. This value of the spin-lattice 
coupling is reasonably small so as to lead to a spin-Peierls transition and not to a structural 
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one 



Now, all the parameters of the Hamiltonian have been determined and we like to check 
the validity of this model by reproducing other experimental results or predicting the value 
of properties still not experimentally measured. Alternatively, the values of J = 160 K and 
A ~ 0.05 could be determined independently by other experiments. 

The first and most obvious check of the consistency of this model is the calculation of 
the spin-Peierls transition temperature Tc. Experimentally, T^. 14 K. To estimate T^ we 
computed the free energies for each lattice size N, and then, we minimized the free energy to 
determine Seq{T, N) as explained above. Tc(N) is the value of the temperature at which the 
dimerization begins. For each lattice size, we used A(N)~^ calculated previously. Results for 
the calculation of 5eg(T,N) for N = 12, 14 and 16 are shown in Fig. These curves have 
some resemblance with the experimental data for the lattice contraction Af, as shown in Fig. 
4 of Ref. 1^. In fact, is related to the atomic displacements u of Eq. (2) and hence to 
6eq- The spin-Peierls transition temperature, calculated for N = 16, is approximately 10.5 K, 
which is reasonably close to the experimental value of 14 K. We have not yet attempted an 
extrapolation of 5eq(T,N) and Tc(N) to the bulk limit. |^ 

As stated above, once 5e(j(T,N) has been computed, any thermodynamical quantity can 
be calculated in the dimerized region. In the first place, we estimated the spin gap as a 
function of T for N = 12, 16 and 20. To simplify this calculation we adopted for the three 
lattices Seq(T, N) corresponding to N = 16. At each temperature, we computed the spin gap 



as the zero temperature singlet-triplet gap of the model with S = Seq(T, N). The results are 
shown in Fig. together with the experimental data from Ref. [Q. Although there are 
strong finite size effects there is a reasonable tendency of the theoretical data towards the 
experimental ones. 

Finally, we have computed the magnetic susceptibility near the spin-Peierls transition, 
for N = 12, 14 and 16. The results are shown in Fig. ^,b. In Fig. it can be seen 
that the susceptibility for each lattice size decays as T goes to zero more rapidly for the 
dimerized Heisenberg model than for the uniform one, which is consistent with a larger spin 
gap in the dimerized case. Notice that for a finite chain there is always a finite gap even 
in the absence of dimerization and for J' = 0. As the lattice size is increased, the spin gap 
in the dimerized model remains finite while the spin gap of the uniform model drops to a 
small value (due to the presence of J'). This behavior explains the fact that the difference 
between the dimerized and the uniform curves is larger as the lattice size is increased, as it 
can be seen in this figure. 

In Fig. we compare the theoretical susceptibility obtained for the 16 site chain with 
experimental data from Ref. obtained with a small magnetic field parallel to the c axis. 
The spin-Peierls transition temperature calculated previously, is shown with an arrow. The 
agreement between theoretical and experimental results is quite good taking into account the 
approximations involved and the strong finite size effects expected in this low temperature 
region. 

In summary, we described the magnetic susceptibility, the temperature dependence of the 
spin gap and the spin-Peierls transition temperature of the CuGeOs using a one-dimensional 
antiferromagnetic Heisenberg model with nearest and next nearest interactions. We obtained 
a quite satisfactory overall agreement with experimental results with only three free param- 
eters. This agreement gives in turn support to the interpretation of the observed features as 
a spin-Peierls transition. The nearest neighbor exchange is approximately equal to 160 K, 
and the ratio of next nearest to nearest neighbor exchange constants is approximately 0.36. 
This value of the parameter a2 would imply a spin gap of the order of 2.4 K, even in the 
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absence of dimerization. If this spin gap is not detected experimentally, it is quite apparent 
that one should necessarily adopt a two-dimensional model. In fact, some small discrepan- 
cies between our model and experiment should be attributable to weak interchain coupling 
as well to a slight spin anisotropy. However, we don't think that it is relevant to include 
these effects before confirmation of the main consequences derived from the present model. 
Details of the calculations and a systematic study of finite size effects are discussed in an 
enlarged version of this report. [^] 

We acknowledge many useful discussions with A. Greco and S. Koval, especially during 
the early stage of this study, and with R. Calvo. 
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FIGURES 

FIG. 1. The magnetic susceptibility of CuGeOs. Experimental curves labelled a, b and c, 
obtained on a single crystal, are from Ref. 111. Experimental curve with dashed line corresponds 



to measurements on a polycrystal. |18| The solid curve is a theoretical one corresponding to the 
Heisenberg model with nearest and next nearest neighbor interactions, with = 0.36, obtained 
numerically on a chain with 16 sites. In the inset, the theoretical susceptibility obtained for lattices 
with even and odd number of spins for T < 30 K are shown. 

FIG. 2. a) Spin gap (in units of J), at T = 0, for the dimerized Heisenberg model with nearest 
and next nearest neighbor interactions as a function of the inverse of the lattice size and for several 
values of b. b) Spin gap in units of J, at T = 0, as a function of b in the bulk limit. The diamond 
indicates the point corresponding to the experimental value of the spin gap. c) Inverse of the 
adimensional spin-lattice coupling constant A as a function of the inverse of the lattice size. The 
dashed line corresponds to a quadratic extrapolation which leads to A"^ 21.3 in the bulk limit. 

FIG. 3. a) bf.q (see text) as a function of the temperature for N = 12, 14 and 16. b) Spin gap as 
a function of the temperature for N = 12, 16 and 20. Experimental data from Ref. Q are indicated 
with solid circles. 

FIG. 4. a) Magnetic susceptibility, in arbitrary units, near the spin-Peierls transition for N= 
12, 14, and 16 as a function of the temperature. The solid (dashed) curves correspond to the 
non-dimerized (dimerized) Heisenberg chains, b) Comparison of the theoretical susceptibility ob- 
tained for the 16 site chain with experimental data from Ref. obtained with a magnetic field 
parallel to the c axis. 
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